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Dead Zone Accretion Flows in Protostellar Disks 
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ABSTRACT 


Planets form inside protostellar disks in a dead zone where the electrical re- 
sistivity of the gas is too high for magnetic forces to drive turbulence. We show 
that much of the dead zone nevertheless is active and flows toward the star while 
smooth, large-scale magnetic fields transfer the orbital angular momentum radi- 
ally outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of 
the dead zone gas to the magnetic fields, despite the rapid recombination of free 
charges on dust grains. Net radial magnetic fields are generated in the magneto- 
rotational turbulence in the electrically conducting top and bottom surface layers 
of the disk, and reach the midplane by Ohmic diffusion. A toroidal component 
to the fields is produced near the midplane by the orbital shear. The process 
is similar to the magnetization of the Solar tachocline. The result is a laminar, 
magnetically-driven accretion flow in the region where the planets form. 


Subject headings: circumstellar matter — solar system: formation — stars: for- 
mation — instabilities — MHD 


1. INTRODUCTION 


The planets of Sun-like stars form from a protostellar disk consisting of 0.01—0.1 Solar 
masses of gas together with 1% dust by mass. The removal of orbital angular momentum 
from the disk is a key to understanding the origins of the planets as it determines the rate at 


which material spirals in to accrete on the star. The magneto-rotational instability or MRI 
Balbus & Hawley |1991)) transfers angular momentum outward and taps the free energy in 


the differential orbital rotation of the disk to drive turbulence and regenerate magnetic fields, 
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at locations where the gas is sufficiently ionized to couple strongly to the fields. However 
the minimum-mass protosolar disk, constructed by adding sufficient hydrogen and helium to 
the planets to give Solar composition on), is mostly too cold for thermal 
ionization. The minimum-mass disk has surface mass density © = 1700(r/AU)~%/? g cm~? 
and temperature T = 280(r/AU)~'/? K at radius r measured in astronomical units (AU), 
and is in vertical hydrostatic balance. Cosmic rays and stellar coronal X-rays are absorbed 
before reaching the midplane, and ionize only the surface layers. Furthermore, rapid recom- 
bination on grain surfaces leads to low abundances of the free electrons otherwise responsible 


for most of the conductivity. Consequently, much of the region where the planets form is not 
subject to the MRI. The region consists of a laminar resistive interior sandwiched between 


two conducting surface layers where magneto-rotational turbulence causes an accretion flow 
Sano et allbood Fromang et al)l2002 Hener & Nelsoxl 2006). The interior 
is not completely quiescent, as waves propagating from the turbulent layers cause hydrody- 
namic stresses and low rates of accretion aa ie Magnetic stresses can 
occur in the interior under favorable conditions, if small grains are removed, for example by 
incorporation into planetesimals, and ionizing cosmic rays reach the disk unimpeded by the 
stellar wind ea bon. However, the resistive interior of the disk retains a weak 
coupling to the magnetic fields even if grains are present and cosmic rays are absent. The 
weak coupling is the topic of this paper. 


2. CRITERION FOR SHEAR-GENERATED MAGNETIC FIELDS 


Magnetic forces can extract angular momentum from the gas if the fields have both 
toroidal and poloidal components. We derive a criterion for the generation of toroidal fields 
by shear acting on radial fields. The Ohmic diffusion quickly eliminates any vertical magnetic 
gradients in the resistive layer, while the conducting layers above and below prevent the fields 
from escaping. The resistive layer can be approximated by an axisymmetric cylindrical flow 
in Keplerian rotation with a uniform resistivity 7. The toroidal component of the induction 
equation reduces to 0B,/dt = —20B, + 7(V?Bg — Bo/r?) where B, and By are the radial 
and toroidal field components in cylindrical coordinates and Q is the orbital frequency. If 
also the toroidal field varies as a power-law with radius, By ~ r”, the field grows through 
the shear term provided r?Q/n > 3(p? — 1)By/B,. We consider the implications for the 
minimum-mass protosolar disk. The magnetic pressure is unlikely to exceed the midplane 
gas pressure, owing to the buoyancy of the fields. The magnetic pressure is a fixed fraction 
of the midplane gas pressure if the power-law index p = —13/8. The resulting criterion for 
shear-generated toroidal fields to grow to about ten times the radial seed field strength, 


v3/(nQ) > 10, (1) 


=3= 


resembles the criterion (Sano & Stone|2002) for magneto-rotational turbulence v4, /(7Q) > 


1. However the shear criterion equation []allows magnetic activity at far greater resistivities, 
since the flattened shape of the disk requires an orbital speed vg much greater than the sound 
speed c,, which in turn is typically greater than the vertical Alfvén speed vaz. The orbital 
speed is a few hundred times the vertical Alfvén speeds found at the midplane in numerical 
calculations of magneto-rotational turbulence TO T. The part of the dead 
zone where the shear criterion holds is better named the “undead zone” because it can be 
reactivated by radial fields siphoned off from the adjacent turbulent layers. In contrast with 
models where angular momentum is removed from the top and bottom surfaces of the disk 


by a magneto-centrifugal wind (Blandford & Payne (1982; (Wardle & Königl(1993), here the 


torques arise from locally-generated magnetic fields. 


3. RESISTIVITY OF THE PROTOSOLAR DISK 


The minimum-mass protosolar disk with well-mixed micron-sized grains has an equilib- 
rium resistivity low enough for magneto-rotational turbulence only in the surface layers, at 
the radii where most of the planets formed (figure [). A deeper layer can generate smooth 
large-scale fields by shear according to equation [] while the midplane within 5 AU of the star 
is almost completely decoupled from the fields. If the grains are removed by incorporation 
into planetary embryos, then turbulence or shear can generate magnetic fields in every part 
of the disk. 


Figure [I] was constructed as follows. The resistivity 7 = 234VT/£e cm? s~! varies in- 
versely with the electron fraction £e = n,./n,, where ne is the electron number density and 
Nn the number density of neutrals OR CTE We compute the electron frac- 
tion by solving a chemical network including ionization by stellar X-rays 
and % Al radionuclide decay, dissociative recombination, charge exchange with metal 
atoms, and grain surface reactions. The recombination is treated using the reduced reac- 
tion network with grains described by (2006). The fraction of the metal 
atoms free to enter the gas phase is chosen to be 190. The resistivity depends only weakly 
on this fraction. The stellar X-ray ionization rate is taken from a fit to Monte Carlo 
radiative transfer results including scattering (igea & Glassgold [1999), scaled to a stellar 
luminosity 2 x 10% erg s-! in 5 keV thermal X-rays to match young Solar-mass stars 
observed in the Orion nebula '2000). The fitted ionization rate is ¢ = 
2.6 x 107° (r/AU)~*(exp[—,/8.0 g cm] + exp[—¥,/8.0 g cm”?)) s?! where ©, and Xp 
are the mass columns lying vertically above and below the point of interest. The fit un- 
derestimates the ionization rate at mass columns < 1 g cm”?, with no significant effect on 


s= 


the locations of the dead and undead zones. The decay of Al in the dust grains yields 
an ionization rate 4 x 10719 s71 (Stepinskill1999). A fraction 0.1% of the radioactive atoms 
is placed in the gas phase, giving a low rate of ionization even with grains absent. We ne- 
glect cosmic ray ionization owing to uncertainty about the extent to which the interstellar 
energetic particles are blocked by the wind from the young star. 


4. MHD CALCULATIONS 


We show results from two 3-D isothermal stratified shearing-box (Hawley et all ligos} 
[Brandenburg et al)1995) MHD calculations of small patches of the minimum- 
mass protosolar disk using the ZEUS code (Stone & Normanl[1999). The first is placed at 
5 AU and includes well-mixed 1 um grains with dust-to-gas mass ratio 1% and a magnetic 
field with a net vertical component of 6 milligauss. The second is placed at 1 AU and has 
no grains and a net vertical field of 30 milligauss. The net magnetic fields are weak, with 
midplane ratios of gas to magnetic pressure 5 x 104 and 3.8 x 10°, respectively. The pres- 
ence of the net fields ensures that the fastest-growing mode of the linear MRI is spatially 
resolved in the calculations. To the net field in both cases is added a part with zero net 
flux that dissipates readily through reconnection. The vertical component is By sin 27x/L, 
the toroidal component By cos2ax/L and the radial component zero, where x is the radial 
position and L the box width. The strength Bo of the sinusoidal field is 0.1 and 1 Gauss, re- 
spectively. The resistivity is calculated using the ionization and chemical network described 
above except that the X-ray ionization rates are directly interpolated in the Monte Carlo 
radiative transfer results and the fit discussed above is not used. The resistivity is fixed at 
its initial equilibrium value at each height in the first run, while the resistivity is allowed 
to vary in space and time through mixing, ionization and recombination in the second run. 
In other respects the calculations are identical to our recent work (Turner et all[2007). In 
particular, the vertical boundaries lie four density scale heights H = c,/Q above and below 
the midplane and are open, allowing magnetic fields to escape leaving the patch of disk with 
a net radial and toroidal magnetic flux. Among the non-ideal terms in the induction equa- 
tion, ambipolar diffusion is important outside 4H. The Ohmic resistivity is greater than the 
Hall resistivity near the midplane at 1 AU while the two are comparable at 5 AU (Wardle 
2007). Under these conditions the Ohmic term has the bigger impact on the amplitude of 
magneto-rotational turbulence (Sano & Stond|2002). Our calculations therefore include only 
the Ohmic diffusion. The domain extends 2H along the radial direction and 8H along the 
toroidal direction, and is divided into 32 radial by 64 toroidal by 128 vertical zones. 


The two calculations give similar results over 130 orbits. A net radial magnetic field 
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Fig. 1.— Location of the dead zone (red) and undead zone (blue) in a cross-section view of 
the minimum-mass protosolar disk. The dead zone gas is decoupled from the magnetic fields, 
while the undead zone is sufficiently ionized for the generation of toroidal fields by shear. The 
remainder of the disk body is well-coupled to the fields and is turbulent through magnetorotational 
instability. The star lies at the origin, the disk midplane falls along the horizontal axis, and solid 
black contours show vertical mass columns of 1, 10 and 100 g cm”?. The pressure in the vertical 
component of the magnetic field is set to 0.1% of the midplane gas pressure. The ionization 
is due to stellar X-rays and the decay of radioactive ?°Al. Recombination on grain surfaces is 
included in the calculation shown at top, where the dead zone extends to 5 AU in the midplane. 
The grains are 1 um in radius and well-mixed in the gas at a 1% mass fraction. With the grains 
removed, there is no dead zone (bottom panel). 


generated in the turbulent surface layers is mixed to the undead zone where it diffuses to the 
midplane. The shear in the undead zone generates toroidal fields (figure 2) of the opposite 
sign, leading to magnetic stresses that transfer angular momentum outward. Lying just 
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inside the undead zone is a dissipation layer about iH thick that receives magnetic fields 
when turbulent motions overshoot the edge of the well-coupled gas. The resistivity is too 
high for magneto-rotational turbulence and dissipates the tangled component of the fields. 
The Ohmic diffusion time H?/ņ at the base of the dissipation layer is about 10 orbits, so 
that the layer acts as a filter delivering to the undead zone a magnetic field that is a recent 
time-average of the large-scale field at the edge of the turbulent layer. The radial field in 
the dissipation layer and at the midplane reverses every 20-50 orbits. The midplane toroidal 
field changes at approximately the rate —30B due to the shear, indicating that local shear 
generation dominates other source terms (Turner et al] 2009. The shearing reverses the 
toroidal field about 10 orbits after the radial field, so that the field lines swing around to 
trail the rotation and the magnetic stress again becomes positive (figure). A related process 
occurs in the Sun, where magnetic fields generated in the convection zone are carried down to 
the tachocline and grow there through differential rotation (Parkerl|1993). The protostellar 
disk differs in that shear occurs throughout the flow, and the turbulence is shut off by the 
resistivity of the gas rather than a transition to radiative energy transport. 


The resistivity at all heights in the disk is typically unchanged as gas mixes between 
the surface layers and interior, because the recombination on the grains is fast enough to 
keep the ionization near its local equilibrium value. However with the grains removed, 
recombination is slow enough for ionized gas to reach the midplane e 
In our calculation with no grains, the undead zone shrinks after 150 years as ionized gas is 
carried down from the surface layers. The flow settles after 200 years into a new steady state 


with an overall accretion rate 20 times greater. The midplane remains magneto-rotationally 
stable (figure B). 


5. EMPIRICAL CONSTRAINTS 


The mass flow rate corresponding to the stresses in the MHD calculations (figure £) is 
sufficient for the disk to accrete on the star within a few million years, the lifetime inferred 
from the disk fractions in young star clusters of different ages ETET The mag- 
netic contribution to the stress is several times greater than the hydrodynamic contribution 
in the dissipation layer, and at times when the fields are strongest in the midplane. The 
time-averaged midplane magnetic and hydrodynamic contributions are comparable. The 
undead zone rather than being inactive has a mean mass flow rate 4% to 61% of that in the 
turbulent layers. The midplane gas moves toward the star on average, but flows away from 


the star immediately after the reversals in the radial magnetic field when the field lines lead 
in the rotation. The root-mean-square midplane field strength found in the calculations is 
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Fig. 2.— Snapshot of the toroidal magnetic field strength at 55 orbits in a resistive MHD 
calculation of a patch of the protosolar disk at 5 AU including well-mixed 1 um grains. The 
undead zone at center is filled with a uniform, 0.1-Gauss shear-generated toroidal magnetic field 
while patchy fields are found in the turbulent layers above and below. The star lies off-page to 
lower left and the disk midplane is horizontal through the image center. 


about 0.1 Gauss at 5 AU and 1 Gauss at 1 AU, in the range 0.1-7 Gauss inferred from the 
relict magnetization of primitive meteorites that were last melted in the asteroid belt during 
the formation of the Solar system (Cisowski & Hood 1991). The agreement is consistent 
with the idea that the planet formation region is a laminar, magnetized accretion flow for 
some part of its history. 


The original dead zone picture with a fixed mass column of accreting material (Gammie 
1996) yields a unique mass accretion rate approximately 1078 Solar masses per year, inde- 
pendent of the stellar mass. In contrast, observations of T Tauri stars show the accretion 
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Fig. 3.— Magnetic accretion stress versus height and time in a resistive MHD calculation of 
a patch of the protosolar disk at 1 AU. The stress is proportional to the radial and toroidal 
components of the magnetic field. It is horizontally-averaged across the patch and is plotted 
using a double-logarithmic color scale, with outward angular momentum transport in red and 
yellow, and inward transport in black and blue. The calculation lasts 300 years and each panel 
shows a 100-year interval. Black curves mark the edges of the undead zone where the resistivity 
suppresses magnetorotational turbulence. Magnetic activity in the undead zone is driven by the 
diffusion of magnetic fields from the turbulent surface layers. After 150 years, the activity grows 
stronger due to an increase in magnetic coupling as ionized gas is mixed from the surface layers 
toward the midplane. 


rate increases with the stellar mass and has a spread of about two decades at Solar mass 
(Hartmann et al. . We suggest based on figure [4] that a similar range of accretion rates 


can be produced by varying the magnetic flux and the abundance of small dust grains. 


The radial transport of material over large distances in the early Solar system is indicated 
by the presence of crystalline silicate grains in Comet 81P/Wild 2 eee or | 
The crystallization requires temperatures of 1000 K that were found near the Sun, while 
much of the comet is icy and formed at temperatures below 200 K in the outer Solar system. 


It remains to be seen whether grains can be carried large radial distances in the turbulent 
surface layers of the protosolar disk without sinking into the undead and dead zones. 
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Fig. 4.— Accretion rate corresponding to the total stress in the two MHD calculations located 
at 5 AU (top) and 1 AU (bottom). The accretion rate per density scale height H is plotted 
against the distance z from the midplane. The results are averaged from 20 to 120 orbits (solid 
black curves at top and bottom) and after mixing reduces the size of the undead zone, from 200 
to 300 orbits (red curve in lower panel). The vertically-integrated accretion rates are 2.2 x 1078, 
6.4 x 107° and 1.5 x 107” Solar masses per year respectively. Shading of the same colors marks 
the undead zone. 


6. DISCUSSION 


We have found that magnetic stresses control the dead zone evolution even with a full 
complement of micron-sized grains, in a minimum-mass protosolar disk ionized only by stellar 
X-rays and radionuclide decay. Consequently, smooth large-scale magnetic fields can occur 
in a wide range of protostellar disks and drive flows in the dead zone before and during 
planet formation as well as afterwards. 
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Models of the early evolution of the solids in protostellar disks depend crucially on the 
dynamics of the gas and will be affected by the presence of a laminar midplane accretion flow. 


The gas drag forces in turbulence cause collisions between the particles, leading to the growth 
or destruction of solid bodies depending on the speeds involved 
2005). In laminar gas, grains settle toward the midplane on timescales much less than the 
disk lifetime ONSI Layered accretion can cause the rapid formation of 
large bodies as particles colliding in the turbulent surface layers grow large enough to fall 
into the laminar interior where destructive collisions are rare (Ciesla 2007). 


Protoplanets migrate away from the locations where they first formed, by exchanging or- 
bital angular momentum with the disk gas through gravitational torques 
197d kada. The rapid sunward migration of the cores of Jupiter and Saturn is a dif- 
ficulty for the standard model of giant planet formation by core accretion (Hubickyj et al. 
'2005). The migration rate depends on the distribution of gas near the planet, which can 
be altered by magnetic forces maa The midplane magnetic pressure in our cal- 
culations reaches a few percent of the gas pressure, sufficient for radial variations to halt 
or reverse the orbital migration. The outer edge of the undead zone is a plausible location 
for radial magnetic gradients large enough to stop the migration, and magnetic effects are a 
potential solution to the problem of the rapid loss of planetary embryos through migration. 


An important question for future studies is how the mass flow rate varies with the 
distance from the star. The fraction of the column coupled to the magnetic fields generally 


increases with the radius. Possibly the inflowing gas will pile up at some location, leading 

to the development of a local radial pressure maximum (Kretke & Li . If so, gas drag 

forces push solid particles toward the maximum, providing favorable conditions there for the 
boa) 


growth of planetesimals 003). The further build-up of material could 
lead to self-gravitational instability and an episode of rapid accretion on the star 
l1996: similar to those associated with intervals of strong mass outflow 
from young stellar objects (Reipurthl|1989). 


Part of this work was carried out at the Jet Propulsion Laboratory, California Institute of 
Technology using the JPL Supercomputing Facility and with support from the JPL Research 
and Technology Development and NASA Solar Systems Origins Programs. 
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